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Abstract 

Our extraction of the pion-nucleon coupling constant from 7rp elastic scattering data 



is outlined. A partial wave analysis (T w < 2100 MeV) is performed simultaneously 
with a fixed-t dispersion relation analysis (T n < 800 MeV). The irNN coupling 
constant g 2 /4ir is searched to find the best fit. The result 13.73±0.01±0.07 (first er- 
ror statistical, second systematic) is found to be insensitive to database changes and 
Coulomb barrier corrections. This value satisfies important elements of low energy 
QCD like the Goldberger-Treiman discrepancy, the Dashen-Weinstein sum rule, and 
chiral perturbation theory predictions of threshold pion photoproduction. 
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I. INTRODUCTION 

In recent times, the pion nucleon coupling constant g 2 jAix has remained an illusive parameter. 
Many groups have worked on the problem using a number of different approaches, yet there is still 
no consensus on its value to the 1% level. This coupling is a very important input parameter in 
low energy QCD and nuclear physics, hence one requires its value to be determined as precisely 
as possible. In one example, the Goldberger-Treiman relation [jl]] connects the nNN coupling 
constant to the well-known weak interaction and hadronic quantities. The deviation in the relation 
has important implications [Q] in low energy QCD. In another example, the Dashen-Weinstein sum 
rule @] connects the coupling constant (via the Goldberger-Treiman discrepancy) to the ratio of 
the strange and light quark masses. Using this relation and the "textbook" value g 2 /Att =14.3, it 
has been argued that the large quark condensate assumption of standard chiral perturbation theory 
may not be valid so that a "generalized" form of the theory may be required [Q]. From only these 
two examples, clearly it is of fundamental importance to pin down g 2 /An . 

Our recent partial wave analyses of txN scattering data up to T w = 2100 MeV (latest publica- 
tion in [|5|]) have included constraints from a simultaneous fixed-t dispersion relation analysis. Our 
most recent (preliminary) solution (SM99 [§]) adds the forward "derivative" E^ dispersion rela- 
tions and the fixed-t C ± dispersion relation to our suite of forward C ± and fixed-t B± ("Hiiper") 
dispersion relations. In these dispersion relations, the coupling constant g 2 is an a priori unknown 
parameter. The nucleon pole (Born) term is a well defined quantity, and extracting the coupling 
constant from the dispersion relation does not involve extrapolations or interpolations. The ben- 
efits of using fixed-t dispersion relations to obtain tcN scattering parameters like the coupling 
constant have been extensively discussed by Hohler in [0j. 

In their influential analysis, Bugg, Carter, and Carter [|]] employed the fixed-t B + dispersion 
relation (see Eqn. ||) and their own partial wave analysis over a narrow energy range (110-280 
MeV) to extract the coupling, obtaining g 2 /An =14.3±0.2. The analysis was not constrained to 
satisfy this dispersion relation, so that value is simply what came out of the data via the partial 
wave analysis. This value was subsequently used in the Karlsruhe analyses, in particular in the 



fixed-t dispersion relation analysis of Pietarinen [[5p used to constrain the partial wave solution 
"KH80" [[Kjp from which the same coupling constant (14.3) was "extracted" using the so-called 
"Hiiper" dispersion relation. To our knowledge, no analyses were performed to test whether other 
values of g 1 jAix gave better results i.e. better satisfied dispersion relations, or better fits to data. 

In the following, we outline our approach to extracting the pion nucleon coupling constant 
from the nN scattering data, where unlike previous analyses involving dispersion relations, it is 
treated as a free parameter to be determined by x 2 minimization. We discuss briefly the method 
in Section 0, followed by the results in Section ||| and some important systematic checks we've 
made in Section p!A| . We close in Section [V| with a summary and some conclusions. 



II. DISPERSION RELATIONS AND THE COUPLING CONSTANT 

The multi-energy partial wave analysis part of our analysis procedure has been described in 
Refs. [ |Tl| ] ■ This part of the analysis has no explicit dependence on the coupling constant g 2 /An 
and so will not be discussed here. The interested reader is referred to those publications. 

The sensitivity to the coupling constant enters through the fixed-t dispersion relations used 
as constraints to the partial wave analysis. We employ the forward (t=0) subtracted C ± (co), and 
derivative E ± {oj) dispersion relations, as well as the unsubtracted fixed-t Hiiper (B±(v,t) ) and 
C + (v, t) dispersion relations. All are implemented from 20 < T n < 800 MeV, and the fixed-t 
relations are also applied over — 0.3 < t < GeV 2 /c 2 . The iso vector constraints can be applied 
to even higher energies, but presently we are having difficulty satisfying the isoscalar C + and E + 
dispersion relations at the higher energies. Since above 2100 MeV we employ the imaginary parts 
from the KH80 analysis [|T0|], and this solution also does not satisfy those isoscalar dispersion 
relations very well around 1 GeV, the incompatability of our solution and KH80 may be at the 
root of this problem. This is still under investigation. The energy range over which the constraints 
limits the covered momentum transfer range to t > —0.3 GeV 2 /c 2 , which is safely within the 
region where the partial wave expansion converges [|7|]. Nonetheless, the kinematic range over 
which the constraints are applied is entirely sufficient for determinations of the ccoupling constant 



since it spans the A resonance region which dominates the dispersion integrals, and where there 
are the most abundant and precise data sets. 

In general, a fixed-i dispersion relation relates the real part of the amplitude at some energy to 
a principal value integral over the imaginary parts at all energies, plus a nucleon pole contribution 
(Born term), and in the case of subtracted relations, an additional subtraction constant which is 
in general energy independent. (For the definitive discussion of dispersion relations, see Ref. [0].) 
This form is demonstrated by the special case of the forward C~ (to) (isovector) dispersion relation 
in Eqn.[l] : 

where //, M are the charged pion and proton masses, to — T^ + // is the total pion energy, 

2 2 

UJb = 2M iS me (unphysical) pion energy at the nucleon pole, Gj^(u) = a" B 3 ^ is the Born term, 
and C~(/i) = 47ro^" + (l + n/M) is the subtraction constant, with the isovector s-wave scattering 
length a 0+ . Refer to Fig. [TJfor a graphical depiction of each term in Eqn. |TJ. 

A. Implementation of Dispersion Relation Constraints. 

The 7riViV coupling constant (appearing in the Born term) and the subtraction constants are 
a priori unknown parameters . Our approach is to treat them as searched parameters to be de- 
termined by minimizing the x 2 goodness-of-fit to the data and the dispersion relations. For each 
dispersion relation the constraint and its x 2 is evaluated as follows. First, the coupling constant and 
all subtraction constants are fixed to some value. With this set fixed, for each iteration in the partial 
wave analysis, the partial waves are used to determine the principal value integral + Born term + 
subtraction constant at a number of equally spaced kinematical points (u> for forward dispersion 
relations, (u, t) for fixed-t) to yield a prediction for the real parts. The real part is then evaluated 
separately using the partial waves. The difference "_Re(from PWA) - i?e(from DR)" is then used to 
correct the real part of the partial waves for the next iteration, and to calculate a x 2 using a desired 
accuracy (presently set such that x 2 /pt. ~ 1 ) as the "uncertainty". This procedure is iterated until 



the solution converges and a minimum overall \ 2 (fit to data + dispersion relations) is achieved. 
This yields the "best" solution corresponding a particular set of dispersion relation parameters. 
The entire analysis is repeated varying these parameters over a multi-dimensional "grid" and the 
minimum x 2 f° r eacn of these solutions recorded. It is observed that e.g. the x 2 versus g 2 /An curve 
is a parabola, and so the parameters for the final solution are determined by fitting quadratics (or 
bi-quadratics to 2-dimensional plots) to these curves and selecting the parameters corresponding 
to the minimum. 

This elaborate procedure has a number of benefits. From it one is able to define a statistical 
uncertainty for g 2 by the variation which changes the overall x 2 by 1 . A systematic uncertainty 
can be estimated from the constancy of g 2 over the applied kinematic range ("extraction error"), 
and from the variation is the \ 2 minimum for the separate contributions of the dispersion relations 
and each of the three charge channels. 

B. g 2 /An from the Hiiper and B + {v t) Dispersion Relations 

When discussing the coupling constant, two dispersion relations merit special consideration. 
One is the isoscalar B + (u t) dispersion relation: 

2v r°° f ImB + {v',ty 



g 2 v\ - v 2 



ReB + (v, t) / dv 

7T Jw\ 



IT Jv\ V 12 — V 2 
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where v = uj + t/AM is the "crossing" energy variable, v B = (t — 2fi 2 )/AM is the energy at 
the nucleon pole, and v\ is the threshold energy. Written in this unsubtracted form, the coupling 
constant is the only unknown parameter. Bugg, Carter, and Carter used this dispersion relation 
[|P and obtained g 2 by inputing in their phase shifts and then "averaging" over a kinematical 
range (Av At) spanning the A resonance. (Refer to Fig. ^|for an example of this method in our 
own analysis). This technique should yield a reliable estimate of g 2 since it is well known [ |T2] ] 
that the dispersion integral is dominated by the first P 33 (A) resonance, and evaluated there it is 
satisfactorily convergent up to a few GeV (below which there still are abundant data). 

The other dispersion relation most useful for determining g 2 is the so-called "Hiiper" dispersion 
relation (see Refs. [0] and [|10|]): 



(v B ± v) 
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where B(0, t) = - J^ 3 dz/ m , ' ^ . It is a clever combination of the invariant £> amplitudes 
such that at fixed t, the left hand side of Eqn. |3| is linear in vb±v with the y intercept independent of 
t and equal to g 2 /M. This dispersion relation was used by Koch and Pietarinen in their influential 
"KH80" analysis paper JTDp . (See Fig. £]for an example from our analysis). Up to a few hundred 
MeV, the dispersion integrals are dominated by the A resonance and suitably convergent up to ~ 
4 GeV when evaluated there. They also have the property that due to crossing symmetry, the "left 
(ys — v) side" is dominated by 7r~p data, and the "right (ub + v) side" by the % + p data. As will 



be shown in Section [III A 2| , one consequence is that extracting g from this dispersion relation is 



relatively insensitive to the Coulomb corrections used in the partial wave analysis. 

HI. RESULTS 

Some g 2 mapping results from our most recent solution ("SM99") [§] are shown in figures 
and [3| . Figure |2| shows the two-dimensional constant x 2 contours of the coupling constant 
versus the subtraction constant in the forward derivative E + (uj) dispersion relation. This contour 
plot was generated using the results of 25 (5x5) solutions spanning the range shown where all 
other dispersion relation parameters were fixed to their optimal values. There is a distinct and 
deep minimum parabolic in shape with negligible correlation between the two parameters. Fitting 
a bi-quadratic to the contours yields g 2 /An = 13.730±0.009, where the uncertainty is statistical 
corresponding the change A% 2 = 1. Very similar results are observed when plotting constant x 2 
contours of g 2 /Aw versus the scattering lengths. 

The statistical uncertainty derived from the x 2 contour mappings is clearly much smaller 
than the overall uncertainty. One way to estimate the systematic uncertainty is to plot the one- 
dimensional curves of x 2 versus g 2 /Air for each of the charge channels, their sum, and the disper- 
sion relation contributions separately, keeping all other parameters fixed to their optimal values. 
This is shown in Fig.^. One sees that all charge channels minimize near 13.73 (13.70, 13.73, 



13.88 for 7r + ,7r~,and CEX respectively) as well as the dispersion relation contribution (13.76). 
This spread gives one indication of the systematic uncertainty. 

Another indication of the systematic uncertainty comes from the values extracted from the 
B + {vt) and Hiiper dispersion relations. The dispersion relations in the solution SM99 were con- 
strained with the g 2 /4:ir =13.73 derived from the x 2 mappings, but fluctuations with respect to v 
and t can arise. The results are shown in figures f| and |5|. One sees an almost negligible energy and 
t-dependence of only about ±0.03 (0.2%) over the full constraint range up to 800 MeV. We refer 
to this as the extraction uncertainty, and see again that it is small with respect to the variations seen 
in Fig.|. 

A. Systematic Checks 

1. Solution with NO Dispersion Relation Constraints 

A number of checks were made in order to gauge other sources of systematic uncertainty. 
One check was to generate a partial wave analysis solution with no dispersion relation constraints 
whatsoever and then use the resulting amplitudes in the dispersion relations to extract g 2 /An. This 
is in effect the approach used by many other prior works, including that of Bugg, Carter, and 
Carter [|]]. This was done to see if the dispersion relation constraints were "pulling" g 2 away from 
a value preferred by the data alone. This turns out not to be the case. Over the same (z/, t) ranged 
used in the constrained analysis, the dispersion relation-free solution yields g 2 /Atx = 13.66±0.18 
(1.3%) from the B + {vt) dispersion relation (see FigJ^), and 13.66±0.07 from the Hiiper disper- 
sion relation (not shown), where only the extraction uncertainties are quoted. Up to 450 MeV, 
where their are ample modern, precise cross section and polarization data from the meson facto- 
ries, these dispersion relations yield about 13.77±0.07. In fact up to ~500 MeV all the fixed-i 
dispersion relations used in the constrained analysis are reasonably well satisfied. This is an im- 
portant observation: to a decent approximation, the low and intermediate energy scattering data 
exhibit the analytic properties expected of them, and so applying dispersion relation constraints to 



a partial wave analysis solution has the effect of "fine tuning" the amplitudes and not drastically 
altering them from their unconstrained state. 

2. Coulomb Corrections 

The dispersion relations must use amplitudes from which all Coulomb contributions have been 
removed (i.e. "hadronic amplitudes"), consequently a systematic contribution to g 2 can enter 
through the Coulomb correction scheme employed. The direct Coulomb and Coulomb phase 



rotation prescription used in our analysis comes from the Nordita analysis [ |14| ] which was used in 
the Karlsruhe-Helsinki KH80 solution [JTQ|]. Our prescription for the Coulomb barrier correction 
has been criticized as being too simple (see e.g. [|T5|]). To test the effect of our Coulomb barrier 
correction scheme on the coupling constant extraction, we made a solution taking the radical step 
of ignoring it altogether. The minimum in the overall \ 2 was found to be 13.70, varying from 
~13.45 for n + p to 13.72 for n~p and 13.67 for charge exchange. Since surely our Coulomb 
barrier correction is more accurate than ignoring it altogether, we conclude that the systematic 
uncertainty due to our Coulomb barrier correction scheme is not significant with respect to the 
other systematic uncertainties. 

It is interesting to note the effect of the Coulomb barrier correction on the Hiiper dispersion 
relation. We reintroduced the correction into the hadronic amplitudes and calculated the dispersion 
relation as before. The result is shown in Fig]7|. The correction suppresses (enhances) the n + p 
(n~p) amplitudes. Since the plot is dominated by n + p (n~p) data on the right (left) hand side, 
the effect is to "rotate" the line around the intercept, leaving it, hence the coupling constant, 
relatively unchanged. One sees in Fig.[7] that the result changes by only 1.3% (13.55 vs. 13.73). 
Clearly this insensitivity to the Coulomb correction makes the Hiiper dispersion relation valuable 
for determining the coupling constant. 



3. Database Changes 

Another source of systematic uncertainty comes from the elastic pion proton scattering 
database. As has been mentioned, it is well known that the A resonance amplitude dominates 
most of the dispersion relations. It was shown in the context of the Goldberger-Miyazawa-Oehme 
sum rule [ JT6J ] that the difference between the Pedroni, et al. [[FT]] and Carter, et al. ]T8|], total 7r ± p 
cross sections results in a 1% change in g 2 , so it is important to study the database contribution 
uncertainties in the context of our full analysis. 

We constructed a solution where the total cross section data of Pedroni, et al., were removed, 
and one where the total cross section data of Carter, et al., the total CEX reaction cross section data 
of Bugg, et al. [|20~1], and the differential n ± p differential cross section data of Bussey, et al. [ p~9| ] 



were removed. We found only a small change of about ±0.04 with respect to our normal solution. 
As opposed to the result in [fl~6|], the effect of a number of fixed-i dispersion relation constraints 
and many data sets (differential, partial total, and polarization) reduces the sensitivity to any single 
measurement 

We constructed yet another solution where all charge exchange data were removed from the 
database. This CEX-less solution results in a best fit coupling of 13.65, which is only 0.6% lower 
than the nominal value 13.73. This is expected, since one sees from FigJJ that the x 2 minimum 
occurs at a larger g 2 /An for CEX than for n^p . 

It should be noted that in the process of analyzing the ix^p scattering data over the years, many 
new data have entered the database, and we have tried many solutions accepting some and deleting 
other data sets. None of these changes have ever caused a large change in g 2 /An , and the result has 
remained stable around 13.73 since at least 1993. We conclude that the database incompatibilities 
that do exist in the current database do not contribute greatly to the coupling constant uncertainty. 

IV. CONCLUSION 

We have outlined our approach to extracting the pion nucleon coupling constant g 2 /An from the 
Tip elastic scattering database using fixed-t dispersion relations. Obtaining the coupling constant 



from these dispersion relations is theoretically unambiguous and relies on the general principles of 
analyticity, unitarity, and crossing symmetry [[7p. The coupling constant was treated as a searched 
parameter to be determined by a least-squares fit to the data and the dispersion relations. From 
our most recent analysis (solution SM99), our result is g 2 /An =13.73±0.01 ± 0.07, where the 
first uncertainty is statistical (corresponding to a change A% 2 = 1) and the second is systematic. 
The latter uncertainty was estimated from the differences in the x 2 minima for each of the charge 
channels and the dispersion relations (Fig]3|), from changes to the irp scattering database, from 
modifications to the Coulomb barrier corrections (Fig.|7|), from a solution with no dispersion rela- 
tion constraints (Fig. ||), and from the variations over the kinematical range where the dispersion 
relations were applied (Figs. |] and |5|). 

Our result is in serious contradiction with the recent results from np backward scattering dif- 
ferential cross section measurements (see e.g. [ [21) ] and other contributions to these proceedings) 
of g 2 /4n ~14.5±0.3, though more in line with the results of the a similar analysis by Arndt, et al. 
[p2|]. The difference is about 9 of our standard deviations. It is very difficult to see how our result 
could be in error by this amount. Nevertheless we are continuing our efforts to update the analysis 
as data new data come in, and to refine our methods in order to extract the most precise value of 
g 2 /4:7T possible. 

It should be noted that the coupling constant result 13.73 resolves the problem with the 
Goldberger-Treiman discrepancy ( [Jl]], [0) which is too large using the larger coupling. Also it has 
recently been shown [ |23| ] that the Dashen-Weinstein sum rule [|3p (corrected for some higher order 
chiral effects) is not consistent with the larger coupling, favouring instead a value near 13.7. It has 
also been shown [ |24| ] that recent measurements of threshold pion photoproduction are more con- 
sistently described in chiral perturbation theory with the smaller coupling than the larger. Along 
with many other recent determinations which have arrived at a lower coupling constant (refer to 
the review p5|], the fact that these important aspects of low energy QCD are more consistently 
described with a coupling constant near 13.7 rather than 14.3 is in our opinion a good argument in 
favour of adopting the lower value as the current standard. 
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FIGURES 
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FIG. 1. The forward isovector C~(lu) dispersion relation broken down into its component parts: the 
Born (nucleon pole) term, the real part of C, the subtraction constant ("Sub. Cons.") and the principle 
value ("PV") integral over a function of the imaginary part of C (refer to Eqn.|l|). The dispersion relation is 
satisfied if the subtraction constant is independent of energy. In our analysis, at several kinematical points 
the squared difference of the real part and the value required to "flatten" the constant is introduced as a \ 2 
penalty function to be minimized. This ensures that once the solution converges, the dispersion relation is 
satisfied to some prescribed accuracy (~ 1%). This procedure is followed for all the dispersion relations 
used in our analysis. 
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FIG. 2. Figure shows contours of constant total \ 2 in the {E t h, g 2 / Air) plane, where E t h is the sub- 
traction constant in the forward (derivative) i? + dispersion relation. The contours were generated from a 
"grid" of 25 solutions where (E t h, g 2 /A-k) was fixed for each. A clear, deep minimum is observed, as it is 
in general for the contour plots generated for all pairs of dispersion relation parameters. 
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FIG. 3. Best-fit x 2 as a function of the coupling constant g 2 /Air, where all other parameters were 
fixed to their optimal (best fit) values. Shown are the total \ 2 ("data+dispersion relation"), and those for 
the dispersion relations, all data, and three charge channels separately. Note that all curves minimize at 
very similar values (~ 13.73), of which the (small) spread is one indication of the systematic uncertainty in 
<7 2 /47r. The bars indicate that A% 2 =1 (statistical) uncertainty. 



16 



0.3 



* 
* 



> 



o.o 



L373 


g2[avg], dg2[rms]= 


13.74 


0.03 






-- 




/ 
/ 


0.0 


-0.0 


-0.0 


0.0 


0.0 


0.0 






o.i 


-0.0 


-0.0 


0.0 


0.0 


0.0 


0.0 


-- 




' 0.0 


-0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


-- 


0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


-- 


/ 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


-0.0 


-0.0 


-- 


0.1 0.0 

/ 


0.0 


0.0 


0.0 


0.0 


-0.0 


-0.0 


-0.0 


-- 


0.0 0.0 


0.0 


0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


0.0 


-0.0 -0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


0.0 


-0.0 -0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 


-0.0 

1 — 


-0.0 -0.0 

1 h- 


-0.0 

1 — 


-0.0 

1 — 


-0.0 

1 — 


-0.0 

h- 


-0.0 

1 — 


-0.0 

1 — 


-0.0 

1 



100 



T lab [MeV] 



600 



FIG. 4. Results for the B + dispersion relation (Eqn. ^) for a solution using the best-fit dispersion rela- 
tion parameters. The figure shows at each kinematical point (T 7rj t) the deviation from the overall average of 
the coupling constant g 2 /Air . The extracted coupling 13.74±0.03 is very uniform. This dispersion relation 
was not one of the constraints used in the analysis. 
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FIG. 5. Result for the "Hiiper" dispersion relation (Eqn. ||) from our best-fit solution. This dispersion 
relation is constructed so that the curves are linear, the y-intercept gives the coupling g 2 /M, and the left 
(right)-hand side of the figure is dominated by it~p {ir + p) data. This dispersion relation was used as one of 
the constraints, so by construction there is virtually no t-dependence in the extracted couplings. 
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FIG. 6. Results for the B + dispersion relation for a solution where NO dispersion relation constraints 
were used. The extracted coupling 13.74±0.08 is very uniform and perfectly consistent with the best fit 
value (see Figj3|). The consistency and uniformity demonstrates that the scattering data by themselves insist 
on the same coupling as the dispersion relations, and have the expected analytic structure. This "free" 
solution satisfies most fixed-t dispersion relations reasonably well (in particular the isovector ones, with the 
sensitive isoscalar relations are less well satisfied), and yield similar results to the constrained solution. This 
implies that the dispersion relation constraints are merely "fine tuning" the partial wave amplitudes and not 
forcing large changes. 



19 



with Coulomb 




13.88 



13.40 



13.55 



-i> 



500. 

400. 

300. 

200. 

100. 

0. 

-100. 

-200. 

-300. 

-400. 

-500. 
-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 

FIG. 7. Hiiper dispersion relation at t = —0.15 GeV 2 /c 2 , but with the Coulomb barrier correction 
reintroduced into the amplitudes. Since the left side is dominated by ir~p data and the right side by ir + p, 
one sees that the amplitudes are enhanced on one side and suppressed on the other, such that the average 
remains relatively constant (13.55 vs. 13.73 = 1.3% change) with respect to the normal solution. This 
"pivoting" around the y-intercept implies that this dispersion relation is insensitive to the Coulomb barrier 
correction. Therefore the difference between the Coulomb barrier correction used in our analysis and a 
more sophisticated form will not cause a significant change in the coupling constant. 
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